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Dynamical fermions by global acceptance steps* 

Francesco Knechtli^ and UUi WolfP 
Institut fiir Physik, Humboldt Univcrsitat, Newtonstr. 15, 12489 Berlin, Germany 

A study of principle is conducted on the inclusion of the fermionic determinant as a Metropolis acceptance 
correction. It is carried out in the 2-D Schwinger model to prepare later applications to the Schrodinger functional. 
A mixed stochastic/determistic acceptance step is found that allows to include some problematic modes in a way 
to avoid the collapse of the acceptance rate due to fluctuations. 
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1. Introduction 

The inclusion of light dynamical fermions in 
a formerly quenched simulation boosts the com- 
putational demand by a very large factor for 
standard Hybrid Monte Carlo (HMC) type al- 
gorithms. This is even true in cases like the 
Schrodinger functional at weak coupling (small 
physical volume), where we expect the fermionic 
determinant to be less influential. We here ex- 
plore the naive idea to propose in such cases gauge 
configurations by an effcient pure gauge update 
scheme and to filter these proposals by an accept 
reject step that leads to the correct full QCD en- 
semble. Although it is advantageous to optimize 
the coupling or even the type of action to be used 
for the proposals, we here focus on the case of 
proposing with just the gauge part of the full ac- 
tion. A more detailed account is given in f. 

We use the two-flavour two-dimensional 
Schwinger model with a non-compact lattice 
gauge field as a test-case with the partition func- 
tion 



D^i{A)\det{Dw + m)\^, 



(1) 



where Dw is the Wilson Dirac operator. The 
gauge action is included in the normalized mea- 
sure 



Sc — Sf 



(2) 



X/J. 



where Sq is the plaquette action and 



(3) 



a gauge fixing term with the backward difference 
operator A*. Periodic boundary conditions are 
assumed and the ^-function eliminates the inte- 
gration over two zero-momentum modes. In Dyy 
the compact gauge variables = exp(jgv4p) and 
with them the coupling strength g enter. In the 
following it is replaced by the dimensionless vari- 
able z = \faL, where a = g'^ /2 is the string 
tension. The gauge part of the action is purely 
Gaussian. Gauge fields distributed according to 
Dfi{A) can be generated in momentum space and 
then Fourier transformed. This amounts to an 
ideal independent sampling or global heatbath al- 
gorithm. In Fign]the unquenched Metropolis ac- 
ceptance rate 



q = 



with 



i J DfiiA) J D^i{A') 
\<lei{Dw + m)\^ w{A, A') 



(4) 



){A,A') ^ min 1, 



\Aci{D'^ + m)[' 
\(iei{Dw + rn)'[' 



(5) 
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is shown for such proposals as a function of the 
mass. The determinants have been evaluated ex- 
actly in this case (data with error bars). The 
critical mass is defined from the spectral gap 
for individual configurations and mc, ct^ are its 
mean and variance with respect to D^[A). For 
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Figure 1. Acceptance with the exact determinant, 
L = 16. 



the range shown no execptional configurations 
are proposed in practice. The circles in the plot 
are obtained by assuming a Gaussian distribution 
for the fermionic part of the action, measuring 
its width and evaluating q analytically for this 
model. Although the cost to compute exact de- 
terminants in higher dimension is prohibitive we 
here see that good acceptance of our global cor- 
rection is possible. In 1 it is shown that q is ac- 
tually an upper bound for the stochastic method 
discussed next. 

2. Stochastic acceptance steps 

An independent new configuration A' proposed 
as a successor to A is now accepted with the prob- 
ability 



wo{A, A'; T]) = min[l, p{Mi])/ p{7])] 



(6) 



where 77 is a pseudofcrmionic random field gen- 
erated with some distribution p (typically Gaus- 
sian) and M = {D[y + m)^^{Dw + rn) is the 
'ratio'-operator. The main cost of this step is 
just an inversion. It leads to a correct algorithm 
due to detailed balance holding on average, 

|detpH^ + m)|2(u;o(A,A';77))„ = 



which is shown by a change of variables in the 
?7-integral. As mentioned before 



— an < 



(8) 



has been derived. It turns out that for most of 
the physical situations shown in Fig^ the value 
of go is smaller than q to such an extent that 
the stochastic algorithm is rendered useless. We 
found that this is related to the spectrum of A/^Af 
whos eigenvalues Xi also occur in the generalized 
problem with {D\y -\- m){Dw + ti)^ on one side 
and its primed companion on the other. Full spec- 
tra are easily computed in our model and typi- 
cally show a structure of two pairs of eigenvalues 
of order g^,g~^ respectively and the bulk of the 
spectrum near one. The occurence of these spe- 
cial modes can be understood in perturbation the- 
ory. By analytically computing {wo{A, A'-, if))^ as 
a function of {A^} one sees that these extremal 
eigenvalues spoil the stochastic acceptance, even 
if J]^ Ai ~ 1 and q is close to one, and then lead to 
go being very much smaller than q. As a remedy 
to this problem we found that detailed balance 
also holds exactly ^ for a partially stochastic 
(PSD) estimation of the determinant ratio with 
the acceptance probability 



(9) 



\dei{D'^ + m)\^w^{A',A;r^))^ 



(7) 



1, W A"^ exp(-77tp(A/tM - i)Prj) 

Here S is the index set corresponding to the four 
extremal eigenvalues Xi (two from either end of 
the spectrum) that caused the trouble. The com- 
plement of the span of the corresponding eigen- 
vectors defines the projector P. Instead of four 
also any other even number of extremal eigen- 
values may be taken symmetrically from both 
ends of the spectrum, and then this family of 
algorithms interpolates between the exact deter- 
minant w and the fully stochastic method with 
Wo. A few extremal generalized eigenvalues and 
-vectors may for instance be computed with Lanc- 
zos or Ritz minimization methods, but this point 
needs further study. 

Results for acceptance rates are found in FigsEl 
and 121 Here a free parameter t has been intro- 
duced, that allows so vary the size of proposed 
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Figure 2. Acceptance vs. stepsize t {z = 2, L 
8, m = rric + 5ac, ^ = 0). 



Figure 3. Acceptance vs. gaugeparameter ^ (z 
1, L = 8, rn — ttLc + 8f7c, t = 1). 



moves A ^ A" = cos(7rt/2)A + sin(7ri/2)A' with 
A' generated by global heatbath. Also these pro- 
posals obey detailed balance with respect to the 
gauge action. Fig|21 demonstrates that for larger 
<-values PSD remains a feasible algorithm. As 
long as the gauge parameter ^ vanishes, simulated 
fields occur in a completely fixed gauge, while for 
growing ^ the proposals are positioned more and 
more randomly on their gauge orbits. Fig|3|shows 
as one may intuitively expect that this freedom of 
gauge mismatch lowers the acceptance rate. Note 
that q with the exact determinant is independent 
of ^. For a locally gauge invariant distribution p 
the Ty-averages of wq and are invariant under 
transforming both A and A' with the same gauge 
function due to the gauge covariance of D\y but 
not under gauge moves of A' relative to A. Hence 
the acceptance rate and autocorrelation times are 
^-dependent although physical results are not. 

3. Conclusions 

We found that at least in the 2-D Schwinger 
model the inclusion of fermions by a global accept 
step is possible. Gauge fields in this superrenor- 
malizablc model are generally smoother than in 
QCD and this is probably important for the suc- 



cess of the method. In sufficiently small physi- 
cal volume with Schrodinger functional boundary 
conditions however a similar situation arises. We 
hope to be able to apply the global acceptance 
method there and maybe find an efhciency supe- 
rior to HMC. This method has already found 4- 
D applications with blocked and hence smoother 
gauge fields entering into the Dirac operator |2I3| . 
Here the use of HMC seems too complicated due 
to the complex dependence of the fermion ac- 
tion on the fundamental (unblocked) gauge fields. 
Also for such applications it seems to be valuable 
to understand the behaviour of the acceptance 
rate in detail and to look for tricks to boost it. 
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